A new cosmic microwave background constraint to primordial gravitational waves 
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Primordial gravitational waves (GWs) with frequencies > 10~^^ Hz contribute to the radiation 
density of the Universe at the time of decoupling of the cosmic microwave background (CMB). 
The effects of this GW background on the CMB and matter power spectra are identical to those 
due to massless neutrinos, unless the initial density-perturbation amplitude for the gravitational- 
wave gas is non-adiabatic, as may occur if such GWs are produced during inflation or some post- 
inflation phase transition. In either case, current observations provide a constraint to the GW 
amplitude that competes with that from big-bang nucleosynthesis (BEN), although it extends to 
much lower frequencies 10~^^ Hz rather than the ~ 10"^" Hz lower limit from BEN): at 95% 
confidence-level, 57gw/i^ < 6.9 x 10~^ for homogeneous (i.e., non-adiabatic) initial conditions. Future 
CMB experiments, like Planck and CMBPol, should allow sensitivities to Qg^/h^ ^ 1.4 x 10~® and 



^gvih^ < 5 X 10 ^, respectively. 
PACS numbers: 98.80.Cq,95.85.Sz,98.70.Vc 

There are many conjectured sources of a primordial 
cosmological gravitational- wave background (CGWB), 
including inflation, pre-big bangtheories, phase tran- 
sitions, or the ekpyrotic model Such backgrounds 
are among the targets of the Laser Interferometric 
Gravitational- Wave Observatory (LIGO), and they will 
be sought with future observatories, such as NASA's 
Laser Interferometer Space Antenna (LISA), the Big 
Bang Observer (BBO), and Japan's Deci-Hertz Interfer- 
ometer Gravitational-wave Observatory (DECIGO). 

The CGWB amplitude is constrained at the lowest ob- 
servable frequencies, ^ 10"^^ — 10~^^ Hz (correspond- 
ing to wavelengths comparable to the cosmological hori- 
zon today) , by large-angle fluctuations in the cosmic mi- 
crowave background (CMB) temperature 2]. Prospects 
for probing lower CGWB amplitudes at these frequencies 
come from future measurements of the polarization of the 
CMB 110]. Apart from a window around IQ-^ - 10"** 
Hz, where the CGWB is constrained by pulsar timing 
[a, [13 , the strongest constraint to the CGWB amplitude 
for frequencies greater than ~ 10~^° Hz comes from big- 
bang nucleosynthesis (BBN) 7]. The lower limit to the 
frequency range is determined by the comoving horizon 
size at the time of BBN. Primordial gravitational waves 
of shorter wavelengths, or larger frequencies, contribute 
to the radiation density at the time of BBN thereby in- 
creasing the expansion rate and thus the light-element 
abundances. Measurements of light-element abundances 
limit the number of additional relativistic species at BBN 
to the equivalent of 1.4 neutrino degrees of freedom 0, 
which translates to a limit to a current CGWB energy 
density n„„h'^ < 7.8 x 10"^. 



The frequency range ~ 10"^^ — 10"^" Hz remains 
largely unconstrained. An upper limit flg^h^ < 0.1 can 
be placed in this frequency range from QSO astrometry 
[11^. It has been proposed that future measurements 
of anisotropy in the global rate of change of observed 



redshifts might someday get down to figw/i ~ 10 ^ [lOj. 

Here we note that recent measurements of the angu- 
lar power spectrum of the cosmic microwave background 
(CMB) that constrain the nonrelativistic-matter density 
ilmh^ to roug hly 10% [13 are, to a first approximation, 
constraints to the radiation energy density at the time of 
CMB decoupling; the constraint corresponds to a limit of 
a few extra neutrino degrees of freedom. From this, we 
infer that the CMB provides a limit to figwft-^ that may 
be competitive with that from BBN, but extends to the 
lower frequencies, ^ 10~^^ Hz, corresponding to wave- 
lengths comparable to the comoving horizon at CMB de- 
coupling |2,1| . This limit therefore improves upon previ- 
ous constraints over the frequency interval 10"^''' — 10"^° 
Hz by four orders of magnitude. 

More precisely, the CGWB behaves as a free-streaming 
gas of massless particles, just like massless neutrinos, and 
therefore affect the growth of density perturbations in 
ways in addition to their effect on the expansion rate at 
decoupling. If the CGWB energy-density perturbations 
are adiabatic (i.e., have the same density distribution as 
other relativistic species), then the effects of the CGWB 
on the CMB/LSS are indistinguishable from those due to 
massless neutrinos. In this case, CMB/LSS constraints 
to the number of massless neutrino species 12] translate 
directly to a constraint to the CGWB energy density. 
If, however, the primordial perturbations to the CGWB 
energy-density perturbations are non-adiabatic, as might 
be expected if they are produced by inflation, pre-big- 
bang models, ekpyrotic, or phase transitions and/or cos- 
mic turbulence (see, e.g., Ref. [iSj), then the CMB/LSS 
effects of the CGWB may differ from those of adiabatic 
massless neutrinos. 

In this paper, we carry out a detailed analysis of 
current constraints to the CGWB amplitude that come 
from current measurements of the CMB power spectrum 
and matter power spectrum. Our calculations of the 
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CMB and matter power spectra include the effects of the 
CGWB on the expansion rate and on the growth of per- 
turbations, for both adiabatic and non-adiabatic initial 
conditions for the CGWB. We include current constraints 
from the CMB, galaxy surveys, and the Lyman-a forest. 
We then forecast how these constraints may be improved 
with future CMB measurements. 
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FIG. 1: Adiabatic: The marginalized (unnormalized) likeli- 
hoods for the CGWB energy density if perturbations to the 
CGWB density are adiabatic. The dotted curve is the result 
obtained using only CMB data. The thick solid curve includes 
galaxies as well as the Lyman-a forest. In all of the aforemen- 
tioned curves, the marginalization is over the nonrelativistic- 
matter density Q.mh'^ , baryon density Qbh^ , scalar spectral 
index Us, power-spectrum amplitude As, the optical depth 
r to the surface of last scatter, and the angle 9 subtended 
by the first acoustic peak (marginalization over 9 essentially 
stands in for marginalization over the Hubble constant). We 
hold the geometry fixed to fiat, the number of neutrinos to 
— 3.04, and the neutrino masses fixed to zero. Finally, the 
dot-dash curve (to the right) shows current constraints from 
the CMB-|-galaxies-|-Lya if we allow for and marginalize over 
nonzero neutrino masses as well. The number of equivalent 
neutrino degrees of freedom (Ng^r) is shown on the bottom 
axis. Homogeneous: same as the left panel, except for homo- 
geneous initial conditions for the CGWB. The arrow indicates 
the 95% CL upper limit Qg^^h^ < 6.9 x 10^^ that we adopt 
as our central result. This is obtained from the analysis that 
includes current CMB-|-galaxy-|-Lya-|-free m^. 

We first consider the case when the CGWB has adi- 
abatic initial conditions. In this case, the effects of the 
CGWB on the expansion history and structure forma- 
tion are identical to those of massless neutrinos. The 
analysis proceeds just as in Ref. We have updated 

this analysis to include new small-scale CMB results, as 
well as constraints from the Lyman-a forest. The CMB 
resuhs we use are from WMAP, ACBAR, CBI, VSA, 



TABLE I: 

CMB Experimental Specifications for Fisher Matrix 
Experiment ^^bcam 

Planck: 7.1 42.2 80.5 0.8 1.4 X IQ-*^ 

5.0 64.8 132.3 

CMBPol: 4.0 1.0 1.4 0.8 5 X lO"'^ 



NOTES.— The beam width, ^beam- (FWHM) is given in 
arcminutes. Weights, (iut,p)~^^'^, a^re in arcminutes ^K. The sky 
fraction is given by /sky The sensitivities to Slgw/i^ are 95% CL 
for homogeneous initial conditions. 

and BOOMERanG, and we use the measurement of the 
galaxy power spectrum from the 2dF Galaxy Redshift 
Survey and the Sloan Digital Sky Survey (SDSS), and 
the Lyman-a forest 14]. We implement the Lyman-a 
constraints following the method discussed in Ref. (isf . 
with minor modifications that were suggested by the au- 
thors. To translate the constraint to the number of extra 
neutrino degrees of freedom to a CGWB energy density, 
we use the relation f2g„/i^ — 5.6 x 10~^, the density con- 
tributed by a single massless-neutrino species. 

Results for adiabatic initial conditions are shown in 
Fig.n A limit at 95% CL of n^^h^ < 3.9 x IQ-^ is ob- 
tained from a combination of current CMB data, galaxy 
power spectrum, and the Lyman-a forest, and under the 
assumption that the number of neutrino degrees of free- 
dom is Ni, — 3.04 and that neutrino masses are free 
to vary. Due to a slight discrepancy between the mat- 
ter power spectrum from the best-fit CMB model and 
that measured in galaxy surveys and Lyman-a forest 
measurements, the addition of galaxy surveys and the 
Lyman-a forest weakens the bound by roughly a factor 
of two. A small CGWB component improves slightly 
the CMB-|-galaxy-t-Lya agreement (cf., the solid curve 
in Fig. although the difference between Ng^ = 
and A^gw = 2 is statistically insignificant. Although not 
shown, we find that the exclusion of the Lyman-a for- 
est weakens the CMB-|-galaxy-f Lya bound only slightly. 
If neutrino masses are assumed to be undetermined, 
then the CMB-fgalaxy-f Lya bound is shifted by approx- 
imately two neutrinos (cf., the dot-dash curve in Fig.[T]l. 
which indicates that there is a degeneracy between the 
neutrino mass and the CGWB. This same trend has been 
observed in Ref. |^^. Note that the bound is improved 
by roughly a factor of 4 if we include only current CMB 
data. 

If the CGWB is produced from quantum fluctuations 
to the spacetime metric during the same superluminal 
expansion that produced primordial density perturba- 
tions (e.g., from inflation, but also from pre-big-bang 
or ekpyrotic scenarios) , then primordial perturbations to 



the CGWB density should be non-adiabatic. In such sce- 
narios, the particle species in the primordial Universe are 
all produced by decay of the inflaton; this is why infla- 
tion produces primordial adiabatic perturbations; i.e., 
the fractional perturbation to the energy density of all 
the particle species are the same. However, gravitational 
waves are produced during inflation by quantum fluctu- 
ations in tensor perturbations to the spacetime metric — 
not through decay of the inflaton. The CGWB should 
therefore not have the same primordial energy-density 
perturbations as the particle species; in fact, in linear 
theory, there should be no primordial perturbation to 
the CGWB amplitude. We therefore re-do our likelihood 
analysis assuming the CGWB has homogeneous initial 
conditions. More precisely, we have chosen to set the 
initial CGWB density perturbation to zero in the con- 
formal Newtonian gauge. With this ansatz, the primor- 
dial curvature perturbation vanishes in the limit that the 
CGWB energy density dominates, as it should; the cur- 
vature perturbation approaches the standard adiabatic 
perturbation in the limit that the CGWB vanishes, also 
as it should. In this limit (appropriate for our analy- 
sis), self-consistency of the perturbation equations de- 
mand that nonzero higher-order moments in the CGWB 
distribution function are induced at early times. Details 
are given in Ref. [l^ . 

If the CGWB is initially homogeneous, then the ini- 
tial conditions for the CGWB perturbations differ from 
those for massless neutrinos. This will affect the growth 
of perturbations, especially at large scales 0|, and the 
degeneracy between the CGWB and massless neutrinos 
is thus broken. The bound to the CGWB then turns out 
to be stronger than in the adiabatic case. Fig. ^ shows 
results for the likelihood for flg^jh^ for different combi- 
nations of current data sets as well as forecasts for the 
likelihoods expected when future CMB experiments are 
included. If the CGWB is produced by some mechanism 
that leaves its primordial density uncorrelated with the 
curvature perturbation — e.g., inflation or perhaps some 
post-inflation phase-transition mechanism — then this is 
the result that should be applicable. We adopt as our 
95% CL upper bound, ^g^h^ < 6.9 x IQ-^, for homo- 
geneous CGWB initial conditions from the combination 
of data from current CMB experiments, galaxy surveys, 
and the Lyman-a forest and under the assumption that 
the number of neutrino degrees of freedom is N^, = 3.04 
and that neutrino masses are free to vary. Note, again, 
that the bound would be roughly twice as strong if we 
were to restrict ourselves only to CMB data. And again, 
although not shown, we find that the exclusion of the 
Lyman-a forest weakens the CMB-t-galaxy-t-Lya bound 
only slightly. 

Our central results are summarized in Fig. |21 which 
shows rigw/i^ versus gravitational- wave frequency. Our 
new constraints are competitive with the BBN constraint 
over the frequency range where both constraints apply. 
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FIG. 2: The gravitational- wave density figw/i^ versus fre- 
quency. The BBN constraint corresponds to a limit of 1.4 
extra neutrino degrees of freedom. We also show our con- 
straints, from current CMB, galaxy, and Lyman-a data, for a 
CGWB with adiabatic primordial perturbations ("Adiabatic 
(current)") and for homogeneous initial conditions ("Homo- 
geneous (current)"), as well as our forecasts for the sensitivi- 
ties if current CMB data are replaced by data from CMBPol. 
Also shown are the reaches of LIGO and LISA. BBO (not 
shown) should go deeper, but primarily at frequencies ~ 1 
Hz. Large-angle CMB fluctuations (also not shown) constrain 
^swh? < 10"", but only at frequencies < 10"^'' Hz. The 
LIGO S3 upper limit is from Ref. 18] and the msec pulsar 
curve is from Refs. 0,01- 



The precise value of the BBN constraint depends on the 
precise constraint to the maximum number of neutrino 
degrees of freedom allowed by BBN. Some authors "2^ 
claim a limit {N^ - 3.04) < 0.2 (at 95% CL), but more 
recent and conservative estimates (that include new ^He 
measurements and the CMB value for the baryon den- 
sity) j^, which we choose to adopt, place the limit at 
{N^ - 3.04) < 1.4, comparable to the CMB/LSS bound 
we have derived. However, our new results apply four 
decades lower in frequency, and provide the strongest 
constraint to the CGWB amplitude over the frequency 
range lO^^^ - 10-^° Hz. 

To forecast the sensitivity of future CMB experiments 
to the CGWB, we have carried out a Fisher analysis 
that shows that when Planck and CMBPol fly, the sen- 
sitivity should be increased by a factor of roughly 10, 
while the BBN constraint may continue to be limited 
by the same astrophysical systematic uncertainties. See 
Table 1 for the experimental specifications used in our 
Fisher analysis. In our Fisher analysis we included the 
improved CMB observations as well as the current galaxy 
and Lyman-a constraints and allowed TOj/ to vary with 
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= 3.04. 

We have not determined precisely the lower end of the 
frequency range for which our bound applies. In or- 
der for the constraint to apply, the gravitational-wave 
wavelength must be within the horizon at roughly the 
time of, or slightly before, recombination. Otherwise the 
waves do not propagate as massless modes. Analytic and 
numerical integrations of the mode equations for gravi- 
tational waves in an expanding Universe (e.g.. Fig. 2 
in Ref. [131 )j indicate that the mode is oscillating when 
kT ~ 10, where k is the wavenumber and r the con- 
formal time evaluated at decoupling. This translates to 
a frequency ~ 5 x 10~^^ Hz. More realistically, the 
gravitational wave will need to oscillate for a while be- 
fore recombination in order to have the effects we have 
considered here. We therefore tentatively estimate 10~^^ 
Hz as the lowest frequency for which our bound applies, 
although the precise value may differ slightly. We leave 
a more precise calculation for future work. 

There is also a slight correction if our bound is applied 
to a scale-invariant spectrum. In this case, the number of 
gravitational-wave modes propagating as massless modes 
changes with time, as more modes enter the horizon. As a 
result, the energy density does not scale with scale factor 
a simply as a"**. This, however, produces only a loga- 
rithmic correction, which is within the theoretical error 
of the treatment we have presented here. 

Finally, we point out that the limit is probably not rel- 
evant for scale-invariant spectra, such as those produced 
by inflation. Those are already constrained to be roughly 
eight orders of magnitude lower in amplitude, at slightly 
lower frequencies ~ 10~^^ Hz, from large-angle fluctua- 
tions in the CMB. However, phase transitions or other 
exotic mechanisms that produce a CGWB at frequencies 
> 10~^^ Hz will now face this new constraint. 
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